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Abstract 



We consider a left-right symmetric SU (4) c x SU (2) ^ x SU (2) R (4-2-2) model with gravity mediated 
supersymmetry breaking. We find that with 4-2-2 compatible non-universal gaugino masses, t — b — r 
j> ■ Yukawa coupling unification is consistent with neutralino dark matter abundance and with constraints 
^ . from collider experiments (except (g — 2) M ). The gluino mass lies close to that of the lightest 
! neutralino, so that the gluino co-annihilation channel plays an important role in determining the 
\ neutralino relic abundance. By relaxing the Yukawa unification constraint we find stau and stop 
CO ■ masses as low as 200 — 220 GeV. We highlight some benchmark points for these cases with 40 < 
§^ '■ tan/3 < 58. 
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1 Introduction 



As a maximal subgroup of Spin(lO) (commonly known as SO (10)), the gauge symmetry SU(4) C x 
SU(2)i x SU(2)r (4-2-2) pQ captures many salient features exhibited by its covering group. Even 
as a stand alone symmetry group, 4-2-2 implements electric charge quantization, albeit in units of 
±e/6, rather than ±e/3 [21 E]- It explains the standard model quantum numbers of the quark and 
lepton families by assigning them in bi-fundamental representations and it also predicts the existence 
of right handed neutrinos pQ. However, there are some important differences between models based 
on 5*0(10) and 4-2-2 which, in principle, can be experimentally tested. For instance, in 4-2-2 the 
lightest magnetic monopole carries two quanta of Dirac magnetic charge [I]. (In SO(10) the lightest 
monopole carries one quantum of Dirac magnetic charge, unless SO(10) breaks via 4-2-2.) By the 
same token, 4-2-2 predicts the existence of SU(3) color singlet states carrying electric charges ±e/2 
[31 [S]. Finally, gauge boson mediated proton decay is a characteristic feature of 50(10) which is 
absent in the 4-2-2 framework. 

While these different experimental signatures can help distinguish 50(10) from 4-2-2, they mostly 
rely on physics operating at superheavy scales. A major motivation for this paper is to highlight 
some important differences in the low energy predictions of supersymmetric 50(10) and 4-2-2 models, 
stemming from the Higgs and sparticle sectors of these models. An exciting new feature is that these 
predictions can be sufficiently different so that they can be compared at the LHC. 

Supplementing 4-2-2 with a discrete left-right (LR) symmetry [H [6] (more precisely C-parity) [7] 
reduces from three to two the number of independent gauge couplings in supersymmetric 4-2-2. In 
combination with Yukawa unification [5], this has important implications for low energy Higgs and 
sparticle spectroscopy which we will explore in this paper, and compare with the corresponding 
predictions from an SO(10) model. 

In 4-2-2 the matter fields are unified into three generations of ip (4, 2, 1), and the antimatter fields 
are in three generations of ip c (4, 1, 2). If the MSSM electroweak doublets come from the bi-doublet 
H(l, 2, 2), the third family Yukawa coupling H ip c ip yields the following relation valid at the GUT 
scale (Mgut), namely 

Yt — Yb — Y T — Y Dirac . (1) 

We will assume that due to C-parity the soft mass 2 terms, induced at Mgut through gravity 
mediated supersymmetry breaking J9|, are equal in magnitude for the scalar squarks and leptons 
of the three families. The asymptotic MSSM gaugino masses, on the other hand, can be non- 
universal from the following consideration. From C-parity, we can expect that the gaugino masses 
at Mqut associated with SU(2)l and SU(2)r are the same. However, the asymptotic SU(A) C and 
consequently SU(3) C gaugino masses can be different. With the hypercharge generator in 4-2-2 given 
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by Y = y|(-B — L) + y |-^3_r, where B — L and are the diagonal generators of SU(4) C and 
SU(2) R , we have the following asymptotic relation between the three MSSM gaugino masses: 

M a = ^M 2 + ^M 3 . (2) 

The supersymmetric 4-2-2 model with C-parity thus has two independent parameters (M 2 , M3) in 
the gaugino sector. 

In this paper we wish to explore whether Yukawa coupling unification in 4-2-2 is compatible 
with recent observations of the dark matter relic abundance and other collider-based experimental 
constraints. A similar analysis for SO (10), which we closely follow, has been carried out by Baer 
et al. [10] . Solutions consistent with SO(10)Yukawa unification have been obtained in [TU] only for 
very special values of the fundamental parameters. Furthermore, it turns out to be quite difficult in 
this model to reconcile the lightest neutralino primordial abundance with the observed dark matter 
densities. 

By introducing non-universality in the gaugino sector, we can allow the neutralinos in 4-2-2 to 
be closely degenerate in mass with the gluino, which is not possible in 50(10). This opens up, in 
particular, the bino-gluino co- annihilation channel [TTJ, which turns out to be an essential difference 
between the 4-2-2 and 50(10) models. In order to make Yukawa coupling unification compatible 
with radiative electroweak symmetry breaking (REWSB), one needs to implement some splitting in 
the Higgs sector, with m 2 Hu < m 2 Hd . Such a splitting may be introduced via a D-term contribution to 
all scalar masses [12], or it can be generated via GUT scale threshold corrections related to a large 
Dirac neutrino Yukawa coupling [13]. It has been noted [H] that a splitting just in the Higgs soft 
terms, as opposed to splitting in all scalar masses, yields better Yukawa unification, and so we focus 
on this approach. Since one of our goals is a comparison of 4-2-2 and 50(10) models, we follow the 
same notation as in [TO]. We parameterize the Higgs soft mass splitting by m% = m\ Q =p 2M|>, 
where m^ is the MSSM universal Higgs soft mass 2 term. The supersymmetric 4-2-2 model we are 
discussing thus has the following fundamental parameters: 

mie,m w , M D , M 2 , M 3 , A , tan f3, sign fi. (3) 

Thus, compared to the 5O(10) model of [10], we have one additional parameter in 4-2-2 which plays 
a crucial role in realizing Yukawa unification consistent with the desired neutralino relic density. 

The outline for the rest of the paper is as follows. In Section [2] we summarize the scanning 
procedure and the experimental constraints that we have employed. We present the results from our 
scan in Section [3j where we compare the 4-2-2 and 5O(10) models and then proceed to highlight 
some of the predictions of the 4-2-2 model. Our conclusions are summarized in Section HI 
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2 Phenomenological constraints and scanning procedure 



We employ ISAJET 7.78 package [15J to perform random scans over the parameter space. In this 
package, the weak scale values of gauge and third generation Yukawa couplings are evolved to Mqut 
via the MSSM renormalization group equations (RGEs) in the DR regularization scheme, where 
Mgut is defined to be the scale at which g\ — g^. We do not enforce an exact unification of the 
strong coupling gz — gi — 92 at Mjut, since a few percent deviation from unification can be assigned 
to unknown GUT-scale threshold corrections [16J. At Mgut, the boundary conditions are imposed 
and all the SSB parameters, along with the gauge and Yukawa couplings, are evolved back to the 
weak scale Mz- The effect of the neutrino Dirac Yukawa coupling in the running of the RGEs has 
been shown in [17] to be significant for coupling values ~ 2. In the 4-2-2 model with t — b — t 
unification, the asymptotic neutrino Dirac Yukawa coupling has the same value as ^(Mqut) which 
is relatively small (~ 0.5). Thus, in the following discussion we will ignore it. 

In the evaluation of Yukawa couplings the SUSY threshold corrections [18] are taken into account 
at the common scale Msusy = ^ m i L m t R ■ The entire parameter set is iteratively run between Mz and 
Mqut using the full 2-loop RGEs until a stable solution is obtained. To better account for leading- 
log corrections, one-loop step-beta functions are adopted for gauge and Yukawa couplings, and the 
SSB parameters rrii are extracted from RGEs at multiple scales TOj = m ; (mj). The RGE-improved 
1-loop effective potential is minimized at an optimized scale Msusy, which effectively accounts for the 
leading 2-loop corrections. Full 1-loop radiative corrections are incorporated for all sparticle masses. 

The requirement of radiative electroweak symmetry breaking (REWSB) [19] puts an important 
theoretical constraint on the parameter space. Another important constraint comes from limits on 
the cosmological abundance of stable charged particles [20]. This excludes regions in the parameter 
space where charged SUSY particles, such as f\ or t\, become the lightest supersymmetric particle 
(LSP). We accept only those solutions for which one of the neutralinos is the LSP. 

We have performed random scans for the following parameter range: 



< 



m 16 



< 20 TeV, 



< 



M 2 



< 1 TeV, 



< 



M 3 



< 1 TeV, 



-3 < A /m w < 0, 

< M D /m 16 < 0.95, 

< m 10 /m 16 < 1.5, 

40 < tan/3 < 58, 



(4) 



with fi > , and m t = 172.6 GeV [21]. 
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We first collected 150,000 points for both the SO (10) and 4-2-2 models. All of these points satisfy 
the requirement of REWSB with the neutralino being the LSP in each case. Furthermore, all of these 
points satisfy the constraint flcDMh 2 < 10. This is done so as to collect more points with a WMAP 
compatible value of cold dark matter relic abundance. Once we identify good regions in paramter 
space, we perform a random scan focused around those regions for the 4-2-2 case. After collecting 
the data, we use the IsaTools package [22] to implement the following phenomenological constraints: 

m-± (chargino mass) > 103.5 GeV [2D] . 

rrih (lightest Higgs mass) > 114.4 GeV [23J, 

m f (stau mass) > 86 GeV [20J, 

m- g (gluino mass) > 220 GeV [20], 

BR(B S -> /i + /0 < 5.8 x 10~ 8 [24], 
2.85 x 10" 4 < BR(b -> s 7 ) < 4.24 x 10~ 4 (2a) [25], 

ttcnuh 2 = O.llllHs? (5<t) [26], 

3.4 x 10~ 10 < Aa M < 55.6 x 10~ 10 (3a) [27]. (5) 

We apply the experimental constraints successively on the data that we acquire from ISAJET. As a 
first step we apply the constraints from BR(B S — > fi + fi~), BR(b —>■ s'j), the WMAP upper bound on 
the relic density of cold dark matter, and the (s)particle mass bounds. We then apply the WMAP 
lower bound on the relic density of dark matter, followed by the constraint on the muon anomalous 
magnetic moment = (g — 2)^/2 at the 3a allowed region. The data is then plotted showing the 
successive application of these constraints. 



3 Results 

Following Baer et at. [10] we introduce a parameter R to quantify Yukawa unification. Namely, R is 
the ratio, 

max(y t , y b , y T ) 
mm(2/ t , y b , y T ) 

so that R = 1 corresponds to perfect unification and a higher value of R signifies a larger deviation 
from unification. 

We next present the results of the random scan. We first compare the 5O(10) model with the 
4-2-2 model in Figures [T] and [2] following the treatment in [TU]. In Figure [T] we plot the results in 
the (R,m w ), (R,tanf3) and (Qh 2 ,R) planes for 5*0(10) (left panel) and 4-2-2 (right panel). All 
of these points satisfy the theoretical requirement of REWSB and correspond to a neutralino LSP. 
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In addition, these points satisfy the various experimental constraints listed earlier. The light blue 
points satisfy the constraints from BR(B S — > yU + /i~), BR(b — ► 57), the Higgs, chargino, gluino and 
stau mass bounds, and the upper bound on the relic density of dark matter from WMAP. Shown in 
dark blue are points that also satisfy the lower bound on Xi dark matter abundance. In Figure [2] we 
similarly present results in the (mio/mxe, R), (Mo/mie, R) and (A /miQ, R) planes for 50(10) (left) 
and 4-2-2 (right). It is quite obvious from the results that, as expected, using just a random scan it 
is quite difficult to realize acceptable Yukawa unification in 5*0(10) consistent with the experimental 
constraints. Ref. [TU] employs a modified scanning algorithm based on Markov Chain Monte Carlo 
(MCMC) to search the parameter space more efficiently. It is shown there that they show that only 
the /i-resonance (light Higgs) channel is available to bring the neutralino dark matter density in the 
right (WMAP) ball park. While this channel does yield acceptable Yukawa unification consistent 
with WMAP, it is more or less ruled out by the lower bound of 114.4 GeV on the the SM Higgs mass. 

In the initial sweep of the 50(10) model around 150,000 points were identified, consistent with 
REWSB and the requirement that is a LSP neutralino. Yukawa unification consistent with the 
experimental data was found to be no better than 40%, even if we ignore the constraint from Aa M . The 
4-2-2 model yields 'good' solutions with Yukawa unification to better than 10%. More concentrated 
searches around such 'good' points have yielded 'near perfect' unification. Such concentrated searches 
were not performed for the 5*0(10) model as they have already been reported in pU] with the 
conclusion that a narrow, almost excluded, light Higgs funnel region is the only one that is viable 
from the point of view of Yukawa unification and dark matter relic density. 

We now focus on the 4-2-2 model, which does much better than the 5O(10) model in terms of 
Yukawa unification and most of the experimental constraints, including the WMAP bounds on dark 
matter abundance. The constraint from (g — 2) M is found to be largely incompatible with Yukawa 
unification (Yukawa unification is worse than 35% if one insists on (g — 2)^). From Figures [1] and [2] 
we find that the following parameter values are preferred: 

m 16 > 7 TeV, 
46 < tan/3 < 48 and 50 < tan/5 < 52, 
0.6 < m 10 /m w < 0.8 and m w /m 16 w 1.1, 
0.3 < M D /m 16 < 0.5, 
A Q « — 2mi 6 and A ps — 2.5m 16 . (7) 

In Table [1] we show a few benchmark points that are consistent with Yukawa unification. Point 
1 displays the spectrum corresponding to essentially perfect unification (R = 1.00). Point 2 gives 
a 'light' gluino (~ 265 GeV) consistent with 'good' unification (~9%). Point 3 has the lightest 
stop (1911 GeV), again consistent with respectable Yukawa unification (~7%). Note that most of 
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the sparticles are rather heavy as a consequence of requiring Yukawa unification. Note that for all 
three benchmark points the lightest neutralino (LSP) relic abundance is compatible with the WMAP 
dark matter bounds. This comes about because of the relatively small mass splitting between the 
neutralino (essentially bino-like) and gluino which leads to efficient co-annihilation [TT] . 

Figure [3] shows plots in the (Af 3 ,m 16 ), (R,M 2 /M 3 ), (M a ,m w /m 16 ), (M 3 , M D /m w ), (M 3 ,tan/3) 
and (M3, Ao/niis) planes for the 4-2-2 model. Color coding is essentially the same as in Figure [IJ 
except that we now also show in red points that are consistent with all experimental constraints 
(except (g — 2) M ) and have Yukawa unification better than 10%. It appears that the points with 
Yukawa unification seem to favor a non-universal gaugino sector, with M 2 > 10M 3 . This ratio is 
higher still if we also require these solutions to satisfy constraints from experiments. This, of course, 
does not mean that solutions with M2 ~ M 3 do not exist, as the latter have been reported in |10j . 
However, this does suggest a statistical preference for solutions with a significant splitting in the 
gaugino sector. 

The PAMELA experiment has reported an excess in the observed positron flux with no corre- 
sponding anti-proton excess [28]. It may be possible to explain this 'excess' in the context of SUSY 
with the lightest neutralino as the dark matter candidate. One explanation invokes a neutralino of 
mass around 300 GeV decaying into positrons via 'tiny' (~ 10~ 13 ) R-parity violating couplings |29j . 
This scenario is consistent with Yukawa unification as we can see in Figure HI 

We have stressed that Yukawa unification seems incompatible with the current experimental 
bound on (g — 2) M . If we do not insist on Yukawa unification in 4-2-2, we can find a much lighter 
MSSM spectrum, which is consistent with all experimental constraints (including (g — 2) M ). This can 
be seen from Figures H] and [5j In Figure H] we show plots in the (m^m^o), (m^m^o), (m^±,m^o) 
and (mg,m^o) planes, with the same color coding as in Figure [3j We also show the unit slope 
line in each plot, thus highlighting the stop co-annihilation region, the stau co-annihilation region, 
and the mixed bino-wino dark matter region. In Figure [5] we show similar plots in the (m^, m^o), 
(mA,m^o), (mg,m^,o) and (mg^m^o) planes. We indicate the A-funnel region with the line rriA = 
2m^o. In Table [2] we present points corresponding to the lightest spectrum found in our investigation 
(disregarding Yukawa unification, but consistent with all experimental constraints). Points 1 through 
5 respectively display the spectrum corresponding to the lightest chargino (133 GeV), CP-odd Higgs 
(284 GeV), gluino (268 GeV), stau (198 GeV) and stop (226 GeV). 

4 Conclusion 

The 4-2-2 gauge symmetry, supplemented by left-right symmetry (C-parity) captures many attractive 
features exhibited by the simplest .SO (10) models. One of these features happens, in some models, to 
be Yukawa unification. We have shown that by relaxing in 4-2-2 the assumption of universal gaugino 
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masses, the resulting MSSM models have rather distinctive mass spectra which can be tested at the 
LHC. Moreover, the primordial abundance of the lightest neutralino in this case is consistent with the 
WMAP dark matter limits, something which is difficult to achieve in 50(10) with t — b — r Yukawa 
unification. We have also studied the implications of relaxing the Yukawa unification condition and 
identified several additional benchmark points which also can be explored at the LHC. Finally, we 
wish to note the recent observation that the little hierarchy problem can be largely resolved in the 
4-2-2 framework [30] • The implication of this for sparticle spectroscopy will be discussed elsewhere. 
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50(10) 4-2-2 
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Figure 1: Plots in the (miQ,R), (tan f3,R) and (Qh 2 ,R) planes for SO(10) (left panels) and 4-2-2 
(right panels). Gray points are consistent with REWSB and Xi LSP. Light blue points satisfy the 
WMAP upper bound on \\ abundance and various constraints from colliders (BR(B S — > 
BR(b — > 57), and (s)particle mass bounds). Dark blue points also satisfy the lower bound on Xi 
density. Green points, additionally, satisfy the constraint from (g — 2) M . 
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Figure 2: Plots in the (mio/wig, R), (M D /mi 6 , R) and (A /mie, R) planes for SO(10) (left) and 
4-2-2 (right). Color coding is the same as in Figure [TJ 
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Table 1: Sparticle and Higgs masses in 4-2-2 model (in units of GeV), with m t = 172.6 GeV and 
/i > 0. Point 1 corresponds to exact Yukawa unification (R = 1.00) while points 2 (3) shows the 
spectrum corresponding to the lightest stop (gluino) with Yukawa unification of 10% or better. Note 
that in each case gluino co-annihilation plays as essential role. 
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Figure 3: Plots in the (M 3 ,m w ), (R,M 2 /M 3 ), (M 3 , m 10 /m 16 ), (M 3 , M D /m 16 ), (M 3 , tan/3) and 
(M 3 ,74 /mi6) planes for 4-2-2. Gray points are consistent with REWSB and Xi LSP. Light blue 
points satisfy the WMAP upper bound on Xi abundance and various constraints from colliders 
(BR(B S — > BR{b — > &y), and (s)particle mass bounds). Dark blue points also satisfy the 

lower bound on Xi primordial abundance. Green points, additionally, satisfy the constraint from 
from (g — 2) ll . Points in red represent a subset of dark blue ones that is consistent with 10% or better 
Yukawa unification. 
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Figure 4: Plots in the (mj, m^o), (mf,m^), (m-±,m^o) and (mj,m^) planes for 4-2-2. Color coding 
same as in Figure [3j Also shown is the unit slope line in each plane. 
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Figure 5: Plots in the (m/^m^o), (m^,m^), (m^m^o) and (m e - L ,m^o) planes for 4-2-2. Color coding 
same as in Figure [31 In the (rriA, Wyo) case we also show the line = 2mjo. 
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Table 2: Sparticle and Higgs masses (in units of GeV), with mt = 172.6 GeV and \i > 0. Points 
1 through 5 correspond to the lightest chargino, CP-odd Higgs, gluino, stau and stop for the 4-2-2 
parameter space given in Eq. (j4l) . Note that these points are not consistent with Yukawa unification 
(R > 2.3), but they satisfy all experimental constraints including the one from (g — 2) M . 
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